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Abstract: We calculate the decay r\ — > 3ir at next-to-next-to- leading order or order p e 
in Chiral Perturbation Theory. The corrections are somewhat larger than was indicated 
by dispersive estimates. We present numerical results for the Dalitz plot parameters, 
the ratio r of the neutral to charged decay and the total decay rate. In addition we 
derive an inequality between the slope parameters of the charged and neutral decay. The 
experimental charged decay rate leads to central values for the isospin breaking quantities 
R = 42.2 and Q = 23.2. 
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1. Introduction 

Since its discovery, the eta particle has been under tense scrutiny, fairly recent reviews are 
[1, 2]. The eta decay to three pions is particularly interesting. It can only happen due to 
isospin breaking. This implies that the decay rate vanishes in the limit of equal of up and 
down quark-masses, ignoring electromagnetic effects. The very first attempts to explain 



the decay [3, 4], considering it to be an electromagnetic decay, resulted in an almost zero 
decay rate, which was in flat disagreement with experiment. 

Later on, a combination of current algebra techniques in SU(3) and the partially con- 
served axial-vector current (PCAC) hypothesis could make a better prediction [5, 6]. This, 
however, underestimated the observed decay rate by a factor of a few. PCAC and cur- 
rent algebra were generalized into Chiral Perturbation Theory (ChPT) and brought into 
a modern form [7, 8, 9]. The lowest order for — > 37r was calculated in [5, 6] and the 
next-to-leading order (NLO) in [10]. The main goal of this paper is to obtain next-to-next- 
to-leading-order (NNLO) expression for this decay. 

We give now a short overview of the present situation, closely following the discussion 



The well-known tree-level result for this decay channel can be derived from current 
algebra or ChPT and has the structure 



The prefactor m u — shows that the decay is isospin violating or S , C/(2)y-violating. The 
magnitude of m u — md determines the size of the isospin symmetry breaking coming from 
the strong interaction itself. A precise determination of this quantity is in principle possible 
here because its size has a direct impact on the decay rate. 

Using Dashen's theorem this factor can be obtained from the physical meson masses 
by removing electromagnetic effects to lowest order. Following the line of the argument 
outlined in Dashen's theorem [12] one arrives at 



Under these circumstances, the theoretical decay width is 66 eV to be compared with 295 eV 
from experiment [13]. One may consider three potential sources for this discrepancy. First, 
the violation of the Dashen's theorem due to the electromagnetic interferences increases 
the value of m^ — m u [14, 15] and therefore the decay width, but this is not sufficient as will 
be discussed more in the discussion section. Secondly, electromagnetic corrections of the 
decay amplitude which are of higher order, order e 2 p 2 , are safely negligible in comparison 
with the strong interactions as pointed out in [16] where the analysis of [3, 4] was brought 



Finally the contribution of the higher order chiral effects must be taken into account. 
Especially since the strong 7T7T rescattering in the <S-wave channel may develop a significant 
correction, see e.g. [17, 18]. The NLO corrections were obtained in [10]. The unitarity 
correction at one-loop level is about half of the total NLO effects. This, of course, confirmed 
the fact that by virtue of a large eta mass, significant rescattering effects in the S- wave can 
occur. The vertex corrections and tree graphs make up the rest of the NLO contribution. 
The coupling constants involved at this level can be rewritten in terms of meson masses 
except for Lg. In [8, 10] Lg was estimated by invoking the OZI (Okubo-Zweig-Iizuka) rule 
and comparing with irir scattering lengths. Their finding for the decay rate of r\ — > 3ir is 



in [11]. 




(1.1) 



B(m d -m u ) = m 2 K o - m 2 K+ - m^ + ra\ + . 



(1.2) 



to NLO. 



160 ± 50 which is still far from the experimental value, however with a large theoretical 
error. 

Given the importance of the unitarity correction at NLO, it was deemed necessary to 
estimate this part of the corrections to higher order. This can be done using dispersive 
methods. In [19], extended Khuri-Treiman equations are used to evaluate the two-pion 
rescattering to the decay r/ — ► tttttt. They achieve a moderate modification, an increase of 
about 14% per cent in the amplitude at the center of the Dalitz plot. Moreover, another 
analysis, based on a somewhat different dispersive method, but also restricting itself to 
two-pion rescattering, represented in [20] suggests also a mild enhancement to the real part 
of the amplitude in the physical region. A more model-dependent analysis of dispersive 
corrections appeared recently [21] relying on combining U (3) x U(3) ChPT and a relativistic 
coupled-channels method, finds agreement with data. 

Given all these, our motivation to perform a full NNLO computation is twofold. First, 
due to a relatively large strange quark-mass, the convergence of three-flavour or SU(3) 
ChPT is an a priori question. The reason hinges on the fact that the ratio M^/M^ is 
much larger than M^/Mp and there are possibly large effects as a result of strange quark 
loops. Three-flavour ChPT is probably less convergent than two-flavour ChPT, see e.g. 
[22] for a discussion. In general, one needs to have several terms available in order to 
check convergence. The situation at present is not fully clear, The results for the vector 
form factors, and 7T7r-scattering have an acceptable convergence, while the results for 
the masses and vr-fT-scattering, seem to converge slower, see [23] and references therein. 
Therefore we would like to check explicitly whether one may treat the strange quark- 
mass perturbatively for this process, namely r] — > 3n. In addition, at NLO the unitarity 
correction provided only half of the total correction. It is therefore also of interest to know 
if the other corrections are important at NNLO as well. This is known to be the case for 
Ki4 by comparing [24] and [25]. 

Our finding shows that the full amplitude up to and including order p 6 corrections 
converges reasonably acceptably but we find larger corrections than in [19, 20]. 

In this paper, we perform the full NNLO calculation of r] — > 3tt in standard three- flavour 
ChPT. We do this to first order in the isospin breaking quantity m u — m^. In addition, we 
perform a first numerical analysis with this expression. We have therefore also estimated 
the order p & coupling constants CJ", using a resonance chiral Lagrangian, assuming that 
vector and scalars mesons saturate the C\. In addition, we derive an inequality between 
the slope parameters. 

The layout of this article is as follows. In Sect. 2 a brief introduction to Chiral Per- 
turbation Theory is provided. Sect. 3 describes how to calculate the 77 — ^ 37r amplitude 
in the presence of mixing, the kinematics for the decay and the form of the amplitude 
at NNLO. We also show the Feynman diagrams and provides references to how we deal 
with the loop integrals and renormalization. A short discussion on our analytic results is 
given in Sect. 4 followed by our estimate of order p 6 low energy constants in Sect. 5. A 
main part of the manuscript is the numerical results and comparison with experiment given 
Sect. 6. We first present a discussion on the amplitude level, Sect. 6.1, then compare with 
the earlier dispersive results, Sect. 6.2. The main comparison with experiment is done in 



Sect. 6.3 for the Dalitz plot distributions and in Sect. 6.4 for the ratio of amplitudes. We 
present the results for the value of R and Q in Sect. 6.5. The App. A contains a discussion 
on the cancellations inherent in a and the derivation of the inequality between the slope 
parameters. The remaining appendices present our NLO expression and the dependence 
on the order p 6 low energy constants. 

2. Chiral Perturbation Theory 

One approach to address Quantum Chromodynamics (QCD) at low energy is the appli- 
cation of effective field theories. In [26] and references therein some basic concepts and a 
few interesting examples can be found for this method. As an effective field theory, ChPT 
is constructed based on the approximate chiral symmetry of the underlying theory (QCD) 
and is an expansion in external momenta and quark-masses, momenta and meson masses 
are generically denoted by p and the expansion is in powers of these. The dynamical degrees 
of freedom i.e. pseudo-Goldstone particles, are manifested as a result of the spontaneous 
chiral symmetry breaking of QCD. Weinberg [7] systematized the use of effective field the- 
ory as an alternative to the current algebra formalism, incorporating naturally the chiral 
logarithms. Gasser and Leutwyler in two elegant papers presented this expansion including 
terms of order p 4 and introduced the external field method [8]. They also formulated the 
extension to three light flavours [9]. They found a substantial correction to the tttt scat- 
tering lengths and effective ranges at this order. Reviews of ChPT at order p are [27, 28]. 
This line of work has been further developed to include p 6 corrections, see the review [23]. 
A more introductory recent review is [29]. 

The full action consists of subterms with a definite number of derivatives or powers of 
quark-masses as shown below 

C ef f = C 2 + C A + Ce + --- . (2.1) 

The subscript indicates the chiral order. Quark masses are counted as order p 2 using the 
lowest order relation m 2 = Bq (m u + m^). The lowest order chiral Lagrangian incorporates 
two parameters and has the form 

£ 2 = ((D^UD^) + ( x rf + U X ] )) 
and the next-to- leading Lagrangian or order p A Lagrangian is given as [9] 

L^D^D^U) 2 + L 2 (D^D u U)(D fl U^D u U) 
+L 3 {D^D fl UD»rfD u U) + U{D^D^U){ X ] U + X U ] ) 
+L b {D^D ll U{ X ] U + U* X )) + L^U + X tf) 2 
+L 7 ( X ^U - X tf) 2 + L s (x*U X *U + X U ] xU ] ) 
-iL 9 (F^ u D^UD u U^ + F^ U D^D U U) 
+L w (tfF*UF L n + H^F^F^ + F L F L n + H 2 { X h) • (2.3) 



(2.2) 



The matrix U £ SU (3) parameterizes the octet of light pseudo-scalar mesons with its 
exponential representation given in terms of mesonic fields matrix as 



U(</>) = exp(iV2<f>/F ) , (2.4) 



where 



V2 VE 



<j)(x) 



n3 + K o 



V2 ' s/6 



2r/ 8 



(2.5) 



K~ K -— , 



The covariant derivative and the field strength tensor are defined as 

DpU = d„U - ir^U + iUlp , F^ v = d„l v - d v l» - i [Ip, Q , (2.6) 

and a similar definition for the right-handed field strength. Here 1^ and represents the 
left-handed and right-handed chiral currents respectively. % ls parameterized in terms of 
scalar (s) and pseudo scalar (p) external densities as x = ^Bq(s + ip). In the process 
discussed in this article we have set s = diag{m u ,md,m s ) and p = 1^ = l v = 0. Finally, 
the notation (A) = Trp (A), the trace over flavours. The SU(3) chiral Lagrangian of order 
p 6 contains 94 operators. For its explicit form we refer to [30]. 

Tree diagrams using only vertices from the Lagrangian £2-, provide the lowest order 
term in the expansion. In general, tree-level diagrams make up the semiclassical part 
of the unitarity of the S-matrix. We thus ought to include loop effects. The infinities 
which arise from one loop diagrams with vertices taken from C2 can not be absorbed 
by renormalizing Fq and Bq or rescaling the fields since these contribute at tree level at 
a different order in p 2 from the one- loop diagrams. Renormalization at higher orders 
contains many subtleties. The procedure used here is discussed extensively in [31, 32]. The 
full divergence structure at order p 4 [8, 9] and p 6 [32] is known. For our calculation, the 
cancellation of the divergences is an important cross-check. The nontrivial predictions of 
ChPT, the so called chiral logarithms, are due to infrared singularities in the chiral limit due 
to the (pseudo)-Goldstone boson intermediate states. Lattice QCD computations support 
this logarithmic behavior when compared with ChPT results evaluated in finite volume 
[33]. 

The part not determined by ChPT, the L\ and higher order LECs, C\ at order p 6 , 
encode the information about higher scale physics which has no dynamical role in our 
effective field theory. A detailed discussion on the determination of the order p 4 LECs is 
given in [9] at order p 4 and at order p 6 in [23, 25, 34]. 

A systematic extension to effective field theories incorporating the resonance fields 
may provide a profound theoretical ground to ultimately underpin the values of the LECs. 
Resonance field methodology takes its original form in Sakurai's hypothesis of vector meson 
dominance. It was worked out at order p 4 by [35]. We will only use it for order p 6 LECs 
in the simplified form discussed in Sect. 5. More systematic approaches at order p 6 exist 
[36] but there are also caveats to be observed from short-distance constraints, both positive 
and negative [37, 38]. 



3. Eta decay amplitude: formalism 




3.1 Matrix-elements in the presence of mixing 

In this section we explain how to calculate matrix-elements by the use of the Lehmann- 
Symanzik-Zimmermann (LSZ) reduction formula in the present of mixing. This is a gen- 
eralization of the discussion in [34] Sect. 2, to the case needed here. 

The scattering amplitude is basically the residue of 
the Fourier transformed Green function in the limit where 
all the outgoing particles go on-shell. 

For the case of rj — ► 7r + 7r~7r , the mixing occurs in 
the two external legs involving neutral particles pion and 
eta as illustrated in Fig. 1. For the decay t] — > 7r 7r°7r° 
mixing is relevant in all four external legs. In [34] two- 
point functions were analyzed in all generality as well as 
amplitudes where one external leg could undergo mixing. 
Here, we reiterate some basic ingredients introduced there 
and will lead to derive a general formula that relates the 
amputated four-point Green function to the scattering 
amplitude order by order in perturbative expansion of 
T] — > 7T + 7r~7r . In the ChPT calculation we only retain 
terms to first order in isospin breaking, we therefore can 
use the relation (3.20) and do not need a more general 
formula for rj — > 7r 7r 7r to all orders in the mixing. 

The amplitude or matrix-element for a process with 
n ingoing and outgoing particles can be expressed in the 
form 



Figure 1: The full four-point 
Green function is represented. 
The Oval stands for the am- 
putated four-point Green func- 
tion and circles indicate the full 
two-point functions. The solid 
lines are external mesons and the 
dashed lines labeled by i and j, in- 
dicate the sum over states implied 
in the two external legs where 
mixing occurs. 



n j™m 



ml) 



(3.1) 



The function Gj 1 ...j n (/ci, . . . ,k n ) is the exact n-point Green function 1 and the coefficients 



Zi are defined as 



Gu(p 2 



m iphys) 



iZi 



(3.2) 



P ~ m iphys 

These are often called field-strength or wave-function renormalization factors. 

We begin by decomposing the full four-point function into the amputated four-point 
function and four full propagators or two-point functions for the external legs. 

Gn(p 2 



G 



1238 



m 



1 phys 



>G 22 (p 2 



m 



2 phys 



m 3ph y s) G 8j(p ~ m 8 p hys )g 12i j . (3.3) 
Subscripts in G1238 are designated for the four external particles 2 namely, n , n 2 , 7r 3 and 



x We consider as operators that are in the Green functions the fields as present in the Lagrangian. The 
formula is valid for all operators with a nonzero coupling to single-particle states. 

2 We use a somewhat generic notation here, 1,2 for the nonmixing external legs and 3,8 for the mixing 
external legs corresponding to ty° and r\. 



7r 8 . A summation over the possible values of i,j is implied in (3.3). Notice that what we 
observe as particles are those in the physical basis. The numbering in (3.3) and the rest 
corresponds to the labeling of the fields in the Lagrangian. We express the amputated 
Green function Qyiij in terms of the contributions at the different chiral orders: 

9m = Q% + Q% + Q% + ■■■■ (3.4) 

In view of the fact that only the two fields associated with the neutral pion and 
eta particles mix, the relevant two-point functions Gij constitute a two-by-two matrix 
G = (Gij) with i,j = 3,8. This can be done since the only nonzero off-diagonal elements 
are G 3 g = Gg 3 . For this matrix we define G _1 = — iP and V can be written as V = P~ 1 +T1. 
P^ 1 = diag(p 2 — rrif ) is the matrix of (inverse) lowest order propagators. We assume here 
that the fields in the Lagrangian have been diagonalized to lowest order. The quantity 
II denotes the sum of all the one-particle-irreducible diagrams as shown pictorially below. 

The function G in matrix form is as follows 

c = g(v 2 ) = * ( Vs8ip2) ~ v ^p 2 A r^i 

[p) r 33 (p 2 )r 8 s(p 2 )-r 2 8 (p 2 ){-r 3 s(p 2 ) r 33 ( P 2 ) ) {6 ' b) 

The denominator is in fact detV(p 2 ) and vanishes at the physical masses given by poles 
in G. One may exploit detV(p 2 ) = to obtain the perturbative corrections to the bare 
mass [34]. By expanding the denominator around the physical mass and knowing that det 
V(p 2 ) vanishes at the limit of physical mass, we can then identify the field-strength factors 
for the neutral pion and eta as 



d 

dp 2 



ml) , 



We define 



(det TV)) 2 2 A (det TV)) 

OIL 



—V 33 (m 2 ). (3.6) 



m iphys 



(3.7) 



dp 2 

and expand also the remaining quantities by chiral order 

i'v • n;r • •••• ^ = i + 4 1) + 4 i) + --- • (3.8) 

Note that by definition II only starts at NLO in the chiral expansion. 

We now have all the ingredients to put in (3.1) and perform the chiral expansion for 
the amplitude. We obtain using Z\\ = Z22 and Tl 3 g,(p 2 ) = Hs 3 {p 2 ) up to order p 6 : 

Ax™ - A {2) 4 A {4) 4 A (e) 4 • • • (3 91 

•/U238 — -^1238 ^ ,/i 1238 ^ • /i 1238 ^ ' \°- J > 

•^1238 = ^1238 1 (3.10) 



rr( 4 ) tt( 4 ) 

am -c (4) - ( z (4) 4. 1 7 {4) 4- V 4) \ c {2) - 38 g {2) _ n 38 c (2) mil 

-'i2.;s i z -ii ^33 1 o 88 1^1238 Am 2 88 Am 2 1233 ' ! - )M ' 



,(6) _ c (6) _ 1 /„ 7 (6) 7 (6) 7 (6)\ c (2) _ 1 / 97 (4) -(4) ~(4)\ c (4) 
-^1238 ~~ ^1238 2 V 11 33 88 / 1238 2 V 11 33 88 / 1238 

+ 8 ((^ 33) ) + ( Z ^) ) ^ 1238 + (^33^) ^1238 

■ ( 1 7 (4) 7 (4) , 1 7 (4) 7 (4) 1 7 (4) 7 (4)\ c (2) 
+ I 7^33 ^88 + 2 11 11 2 11 11 / 1238 

n 3 8(3) (4) ^ ( 4) n 38 (8)( 4 ) ^( 4 ) n 38 (3)( 6 ) ^ (2 ) n 38 (8)( 6 ) „ (2) 

Am? ^ 1288 Ami 1233 ^? 1288 Am 2 ^ 1233 



°2 

n 38 (3)W n 88 (3)W ^ (2) n 38 (8)Wn 88 (3)( 4 ) ^ (2) 

Am 2 ^ 1288 + Am 2 ^ 1233 




7 ( 4) n 38 (3)( 4 ) v (4) n 38 (3)W (4)1138(3)^ \ (2) 

J33 Am? +Z§8 Am 2 + ^ Am 2 J y 1288 



( 4) n 38 (8)W ( 4 ) n 38 (8)W (4) n 38 (8)W \ (2) 

^ 33 Am 2 + Z§8 Am 2 + Al Am 2 J &1233 ' (3 ' 12) 

We defined the abbreviations Am 2 = m 2 — m 2 and Am 2 , = m 2 — m 2 with first a physical 
mass and the second term the lowest order mass. are evaluated at the physical 

7T° mass for I = 3 and at the physical i] mass for 1 = 8. The are evaluated at the 
physical charged pion mass for the legs with indices 1 or 2 and at the physical 7r° mass 
for the leg labeled i and at the physical rj mass at the leg labeled j. The set of terms are 
shown at the different orders. 

3.2 Kinematics and isospin 

We write the amplitude for the decay r)(p v ) — > 7r + (p + )-7r~ (p_)7r°(po) using the Mandelstam 
variables 

s = (p+ + P-f = (P v ~ Po) 2 , 
t = (P+ + Po) 2 = (p v ~ P-f , 

u = (p_ +p ) 2 = (p v -P+) 2 , ■ (3.13) 

These are linearly dependent 

s + t + u = m 2 + m 2 _ + m 2 + + m 2 = 3s . (3.14) 

Due to the SU(2)y symmetry breaking the isospin basis and physical basis in the 7r°- 
77 subset do not coincide. To diagonalize the mass matrix and consequently the two- 
point functions at leading order we perform the following transformation and find the 
corresponding lowest order mixing angle 

7r 3 = 7rcos(e) — ?7sin(e) 

i]s = 7rsin(e) + ??cos(e) (3.15) 

The lowest order mixing angle is 

\/3 m d - m u 



tan(2e 



2 m s — rh 

rh = (m u + m d )/2 . (3.16) 



G-parity requires the amplitude to vanish at the limit m u = and therefore it must 
inevitably be accompanied by an overall factor of m u — which we have chosen to be in 
the form of sin(e). 

A(rj -► tt+tt-tt ) = sin(e) M(s, t, u) (3.17) 
Since the amplitude is invariant under charge conjugation we have 

M(s,t,u) = M(s,u,t). (3.18) 

Note that the isospin breaking factor which is pulled out is different in different references, 
[10] and [39] use different quantities. We have chosen sin(e) since it is the factor that 
naturally shows up at lowest order. 

For the eta decay to three neutral pions, the amplitude must be symmetric under the 
exchange of pions (Bose symmetry) and this together with isospin symmetry and using the 
fact that the decay is caused by the AI = 1 operator (1/2) (m u — m,d) (uu — dd) implies 

A(rj -► ttW) = sin(e) M(s, t, u) (3.19) 

with 

M(s, t, u) = M(s, t, u) + M(t, u, s) + M(u, s, t) . (3.20) 

This relation is only true when isospin breaking in the amplitude is taken into account to 
first order only. M(s,t,u) and M(s,t,u) are treated in the isospin limit. We will work in 
this limit in the remainder of the paper. 

3.3 A simplified form for the amplitude to order p 6 

The scattering amplitude can be represented in terms of components with definite isospin 
assignments as [20] 

M(s,t,u) = M (s) + (s - u)Mi(t) + {s- t)-Mi(u) + M 2 (t) + M 2 (u) - \m 2 {s) . (3.21) 

The function Mj(x = s,t,u) indicates scattering in the kinematic x-channel with total 
isospin I. Analyticity, unitarity and crossing symmetry as imposed on the S-matrix, give 
rise to this exceptionally useful representation. For the derivation and detailed discussion 
for the case of inr we refer to [40]. This relation only holds up to 0(p 8 ). The argument 
is based on the fact that nonpolynomial dependence on s, t, u is related to an absorptive 
part via unitarity. Up to order p 6 there are only absorptive parts in the two-pion S and 
P-waves and then using isospin one derives the form (3.21). 

It is important to note that the polynomial part of the amplitude cannot be split 
uniquely into M/ functions since the relation s + 1 + u = 3so allows a different redistribution 
of the said part to the Mjs. A list of the equivalent redefinitions for the case of K — > 3ir 
can be found in App. A of [41]. The choice we made is to remove as much as possible first 
out of M 2 and then out of Mi. 

Eq. (3.21) makes the formidable task of handling two loop expressions much more 
manageable and we indeed confirmed explicitly the validity of Eq. (3.21) for the decay 
7] 7r°7r+7r- at 0{p 6 ). 
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Figure 2: The Feynman diagrams of order p 2 , (a) and of order p 4 , (b)-(d). 
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Figure 3: Feynman diagrams of order p 6 . 

Note that the neutral decay amplitude can also be expressed directly in terms of the 
Mj(t), 

M(s, t, u) = M (s) + Af (t) + M (u) + \ (M 2 (s) + M 2 (t) + M 2 (n)) . (3.22) 
when using the isospin relation (3.20). 

3.4 Feynman Graphs 

We have collected all the amputated Feynman diagrams needed. In these figures a filled 
circle denotes a vertex of order p 2 , a filled square a vertex of order p and the grey filled 
square a vertex of order p e . 

3.5 Regularization, renormalization and integrals 



As a regularization method we use dimensional integration. The regularization method is 
described in detail in [31, 32]. 



The one-loop functions needed are denned in many places. We use the definitions as 
given in [42, 25, 43]. The two- loop integrals we evaluate numerically using the methods 
described in [42] for the sunsetintegrals and in [43] for the vertexintegrals. 

In addition to the integrals given there, we also need derivatives w.r.t. one of the 
masses in order to be able to pull out the overall isospin breaking factor in the amplitudes. 
For all the one-loop integrals needed, these derivatives can be derived exactly. An example 
of a relation we derived, needed to obtain agreement with the order p 4 result of [10] is 

7=7/ 2 2 2 2\ @ ~Dt 2 2 2\ 

C(m ,m ,m ,p ) = — — ^B(m 1 ,m 2 ,p ) 

dm l m 2 =m 2 =m 2 

= - p T^{B{rn\ m \p 2 )-B( m \ m \Q)-^ . (3.23) 
For the two-loop integrals, we have taken the derivatives numerically. 

4. Analytical results 

4.1 Order p 2 

At leading order there is one tree graph from £2 to compute since we have already diago- 
nalized the fields in the lowest order Lagrangian. The lowest order decay amplitude takes 
on the form 

M {2 \s) = ± (\ml - s) (4.1) 

using the identity s + t + u = 3so- and m n are the physical pion decay constant and the 
physical pion mass respectively, which involve corrections of order p 4 and p G . The higher 
order corrections due to the redefinition of the parameters are carried to the respective 
amplitude of order p 4 and p 6 . This agrees with the known lowest order results [5, 6, 10] 
but is written in a somewhat different form. We have chosen this form since it brings out 
the Adler zero explicitly. 

4.2 Order p 4 

At this order we obtain vertices from £2 to construct the tadpole, Fig. 2(c) and the so-called 
unitarity graphs Fig. 2(d) and in addition, the tree diagram from £4, Fig. 2(b). These are 
the diagrams contributing to G[%8 of (3.11). The second set of terms in -4^38 m (3-H) is 
from what is usually called wave function renormalization. The final two terms in ^4^238 °f 
(3.9) are what are called the mixing corrections 3 . 

The sum of all the contributions described in Eq. (3.11) gives the full one- loop am- 
plitude. G-parity requires the amplitude to be proportional to sin(e) as was the case at 
leading order. This does not occur explicitly at order p 4 , because now isospin is broken in 
the meson loops. In order to pull out an overall mixing angle we then carry out a Taylor 



3 In our calculations, lowest order mixing is dealt with exactly. The mixing at higher orders is dealt with 
perturbatively. 



expansion of the loop integrals involving two charged kaons around the neutral kaon mass. 
The amplitude is written as 



M^(s,t,u) = JL 



M,S 4) (s) + (s - u)M[ A) {t) + (s- t)M[ A \u) + M 2 (4) (t) + M 2 (4) (u) 

(4.2) 



The full expression at 0(p 4 ) for the M (4 \t) can be found in App. B. There, all the masses 
and the pion decay constant are the physical ones. We have taken the expressions at lowest 
order and one-loop order and added the correction terms needed at the higher orders to 
bring them into the form we show in (4.1) and App. B. and are corrected up to order p 6 . 
We retain all the p 4 effective constants in our expression, in contrast to [10] where all but 
L 3 are eliminated in favor of measurable quantities. We have also found our analytical 
results in full accord with that in [10]. 

4.3 Order p 6 

We again split the amplitude as in (3.21). 



M^\s) + {s- u)Mf\t) + {s- t)Mf\u) + M 2 (6) (i) + Mf\u) 

(4.3) 



However, the order p 6 expression is extremely long. The dependence on the order p 6 LECs 
is given in App. C. We split the result in several parts 

Mf\t) = Aff(t) + M\ L (t) + Mj{t) . (4.4) 

Mf{t) contains the contributions from the order p 6 LECs, M^ L {t) contains the contribu- 
tions that involve the order p 4 LECs and Mj(t) is the pure two- loop contribution, only 
dependent on the masses of the pseudoscalars. Mj{t) itself we split in the parts coming 
from vertex- integrals, sunset-integrals and the rest. The latter split is definitely not unique. 
It depends on how many of the relations between the various integrals are actually used 
and we observe strong numerical cancellations between its different parts. We therefore 
only quote numerical results for Mj(t) whole. 

The calculation has been performed independently by each of the authors, the diver- 
gences agree with those of the general calculation [32] and nonlocal divergences cancel as 
required. We have also checked explicitly that the amplitudes can be brought into the form 
(3.21). The numerical results are also done twice and we have checked that they agree. 
Finally, //-independence has been checked numerically and found to be satisfactory. The 
latter is not exact, since we have expressed the order p 4 and p e in the physical masses and 
there is a residual p s \i dependence left over. We do observe a strong cancellation between 
the order p 4 and p 6 /i-dependence as expected. 



5. Resonance estimates of the p 6 LECs and the inclusion of the rf 

Chiral symmetry imposes no constraints on the values of the low energy constants, nev- 
ertheless these constants do depend on the parameters of the underlying theory, QCD, 
namely the masses of the heavy quarks and the scale Aqqd- Hence, all the LECs may 
be determined from first principles employing the Lattice QCD technique. At order p 4 , 
most of the LECs have been determined phenomenologically and some of them are checked 
numerically by the use of the Lattice QCD [44]. At the present time, only very few of 
the order p 6 LECs are estimated using available data[23]. In this section, we will dis- 
cuss briefly how to do an approximate estimation of the LECs within the framework of 
the resonance effective field theories. In the limit N c — > oo, the matching of the QCD and 
ChPT might become feasible since now, an infinite tower of massive narrow hadronic states 
emerge. With these states one can construct a chiral invariant Lagrangian incorporating 
both Goldstone mesons and resonance fields. In practice one has to do a truncation on 
the hadronic spectrum and limit the resonance multiplets to low-lying excitations. The 
construction of the p 4 resonance Lagrangian is discussed in the pioneering works [35, 37]. 
The extention of the earlier works to 0(p 6 ) is presented in [36]. In the following we present 
a resonance Lagrangian at order p 6 with the vector realization of the vector fields, [35, 36] 
use the antisymmetric tensor formulation. The exchange of axialvector resonances does 
not contribute to the process i] — > Sir and they are not discussed here. We do, in addition 
to the lowest vector nonet, include the pseudoscalar singlet and a scalar nonet. 

The building blocks one needs for the construction of the Lagrangian take the form 

u » = iu^D^Uu^ with u 2 = U, 

X± = U ^X^ ±ux*u. (5.1) 

These transform as an octet under SU(3)v using the general CCWZ construction 4 . 

The resonance fields are counted as order 1 in the chiral counting. We do not include 
the complete possible list of terms here but restrict to a smaller subset which contains the 
lowest order interactions of the resonances in our chosen representation for them. This is 
the same subset of possible terms used in most of the work on NNLO ChPT. We only show 
the terms relevant for 77 — > 37r here. For the vector fields we use as Lagrangian 

Cy = -\{V, v Vn + \m 2 v (V^) - -^=<V>KX]> + / x (^Kx-]> • (5-2) 

V/xi/ = V M V^ — V U V^. The singlet component does not contribute to order p 6 for r\ — »■ 3ir. 
For the scalar meson nonet, the matrix of fields S, we consider 

C s = \(V^SV^S - M 2 S S 2 ) + c d (Su^) + c m (S X+ ) (5.3) 

At tree level integrating out the heavy fields is equivalent to applying the equation of 
motion for the elimination of the heavy fields. To solve the equation of motion we perform 



4 See Refs. [35, 36] for a more extensive discussion and references 



a perturbative expansion of the resonance field with coefficients in increasing powers of 
1/Mr {Mr, the resonance mass) and then solve the equation of motion recursively. We 
obtain for the Vector V 1 and Scalar S resonance field 

s = W/"^ + + ^W') + ^(*^*+> + - ■ (5-4) 

Substitution of (5.4) in the Lagrangians (5.2) and (5.3) we obtain the effective action from 
V,S-exchange[31, 25] 

Cv = -i^((^ V ^K,^]-/x^K,X-]) 2 ) (5.5) 

Cs= 2M^ (^ V "M + cm^x+f) (5.6) 
s 

The resonance couplings were determined in [31, 25, 35]. The values we use are 

f x = -0.025, g v = 0.09, c m = 42 MeV, c d = 32 MeV, (5.7) 
and for the masses we use 

m v = m p = 0.77 GeV, m s = 0.98 GeV. (5.8) 

The rj plays a significant role in processes involving r\ due to the r\ — rf mixing. Within 
the quark model, rf and i] mix because of the SU(3) symmetry breaking. In the chiral 
limit the pseudoscalar octet becomes massless but the rf has a residual mass as a result 
of the axial 17(1) anomaly. In the combined chiral and large -/V c limit, however, a nonet of 
Goldstone particles emerge and this provides a perturbative framework to investigate the 
dynamical interplay between rf and Goldstone particles, see e.g. [45, 46]. 

One important result is the saturation of the low energy constant L-j by the rf exchange[9] 



d 2 



L V> = _^rn , g 

7 2Ml 



We therefore perceive the rf dynamical effects at order p 4 through its contribution to 
the effective constant Lj. In the light of this realization, the 77 — 77' mixing effect on the 
Ci involved in the decay 77 — > 37r can be obtained by constructing an appropriate U(3) 
Lagrangian at order p 6 . That 77 — 77' mixing can be treated perturbatively for the decay 
rj —>■ 37T is discussed extensively in [45]. We take for the singlet degree of freedom Pi the 
simple Lagrangian 

Crf = l ~d ll P 1 d»P l - \Ml,Pl + idmP^X-) . (5.10) 

Integrating out Pi leads to the order p 4 term with L7 of (5.9) and the order p 6 Lagrangian 
A/ = -^^(x-)^(x-) with d m = 20 MeV. (5.11) 
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Table 1: The resonance exchange estimates of Ci contributing to rj — > Sir. The vector exchange 
results are taken from [48] scalar exchange from [36]. For the singlet eta contribution, see text. 
The result for the C\ is the top row multiplied by the coefficients given in the table. Only nonzero 
coefficients that also contribute to 77 — > 37r are given. We use here the dimensionless version of the 
C\. Only terms relevant for r\ — > 37r are shown. 



The latter can be rewritten in general in terms of the basis of operators of [30]. The result 



is 5 



^(x-)^(x-) = Ois + \o 19 - io 20 + ^0 21 + 20 27 + ^0 3 i - lo 32 + ^0 33 

-2O35 + O37 - ^0 94 . (5.12) 



We have derived the contribution to r] — ► 3n in different ways. An option is to evaluate 
explicitly the resonance exchange directly from the Lagrangians with resonances. The 
second method is as described above, to evaluate the exchange in general in terms of an 
effective Lagrangian of the pseudoscalars only and then calculate with that one. The third 
option is to rewrite the contributions from vector and resonance exchange in terms of the 
order p 6 ChPT Lagrangian [32]. The latter can be done for the scalar octet using the 
results of [36], for the vector nonet using [48] and the pseudoscalar singlet using (5.12). We 

5 This was derived by the authors of [25] but not included in the final manuscript. It agrees with the 
expression shown by Kaiser [47]. 
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Figure 4: (a) The order p 2 , p 4 and p 6 contribution to — ReMo(s). (b) The order p 6 contribution 
to — Mq{£) split into its parts, the contribution from vertex-integrals (V), sunsetintegrals (H) and 
the remaining pure two-loop part as well as the independent part. 



have checked that the first method also gives the Lagrangians (5.5) and (5.6) and that the 
2nd and third method in the end give the same contribution to r\ — > 3tt. The results for 
the Cj are quoted in Tab. 1. In the numerical estimates we have used Fq = F n . 

6. Numerical results 
6.1 A first look 

In this section we present some plots of the amplitude to give a first impression of how the 
higher orders look like. We start with figures showing different contributions to 2(2). 
It should be noted that since terms can be moved around between the Mj(s) these figures 
cannot be used to check whether we have convergence or not. They are shown for illus- 
tration only. We have actually plotted the negative of the quantities defined earlier, this 
makes the size of the corrections easier to judge. 

The input values we use are the physical eta mass, = 547.3 MeV, the average 
Kaon mass removing electromagnetic effects, mx = 494.53 MeV and a pion mass such 
that s + t + u = m 2 + 3m 2 is satisfied for the charged and neutral decay. So we use 
3m^ = 2m^ + + m^ for the charged decay and m 2 = m^ for the neutral decay. More 
general plots were always done with the mass as for the charged decay. The order p 4 LECs 
are set to the values for fit 10 of [34] and we have set the order p % LECs equal to zero. The 
subtraction scale is fi = 770 MeV. 

In Fig. 4(a) we show the contributions of order p 2 , p 4 and p® to Mo(t). One sees an 
acceptable convergence in the physical domain for the decay. In Fig. 4(b) we show how 
the different parts contribute. As one can see, there are sizable cancellations in the order 
p 6 contribution. In Fig. 5(a) we show the various contribution to the absorptive part of 
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Figure 5: (a) The order p 2 , p 4 and p 6 contribution to — ImMo(s). (b) The order p 4 and p 6 
contribution to — ReMi(s). 
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Figure 6: The amplitude AI(s,t,u) along the line t = u. The vertical lines indicate the physical 
region, (a) Shown are the real and imaginary parts with all parts summed up to the given order, 
(b) Same plot but the contribution from the CI has been removed. 



Mo(s). We see that the total p 6 is about the same size as the order p 4 one. To be noted 
is the three particle cut that contributes first at order p e . This allowed Mo(s) to have an 
absorptive part also below the two-pion threshold. This cut gets contributions from the 
vertex- and the sunset-integrals, diagrams in Fig. 3(f) and (g). As expected, this cut gives 
a rather small contribution. We also show a case where the convergence looks bad, M\(s) 
shown in Fig. 5(b). When we look at the full amplitude this large correction is not visible. 
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Figure 7: The amplitude M(s,t,u) along the line t = u. The vertical lines indicate the physical 
region, (a) Shown are the real and imaginary parts with all parts summed up to the given order but 
the contribution from the L\ and C\ have been removed, (b) Same plot but showing the various 
parts. 



It seems to be tempered sufficiently when all the different M/(s) are summed as in (3.21) 
To see this, let us look at the ampli- 

1 L, fit 1 and C, 

tude along the line u = t as a function 
of s for the full amplitude. Here we plot 
first our full result. In Fig. 6 we show 
the order p 2 , the sum of order p 2 and p 4 
and finally the sum of order p 2 , p 4 and 
order p®. Note the shift in the Adler 
zero in going from order p 2 to order p 4 . 

As a final part here we show the de- 
pendence on the subtraction constant [i. 
As mentioned above, our full result is 
/U-independent to the order it should be 
with a cancellation between the varia- 
tion at NLO and NNLO. The ^-depen- 
dence creeps in via the estimate of the 
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C\ and at which scale fi it is applied. In 
Fig. 8 we plot the real part of the am- 
plitude to NLO and NNLO at fj, = 0.6, 
0.77 and 0.9 GeV. The variation with \i 
is fairly small. 
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Figure 8: The amplitude M(s,t,u) along the line 
t = u to NLO and NNLO order for three choices of 
the subtraction point /x, i.e. three choices at which 
the estimate of the Cl(fi) is applied. 
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Figure 9: (a) Decay amplitude obtained by use of extended Khuri-Treiman equations [19] along 
the line s = u. (b) Alternative dispersive analysis for the decay amplitude [20]. Figs, from [49], 
adapted from [19, 20]. 



6.2 Comparison with the dispersive result 

Since long ago it has been known that the decay amplitude of rj —* iririr receives a sizable 
enhancement due to the loop correction at 0(p 4 ) [10]. In fact it turned out that a large 
part of this correction comes from irir rescattering in the final state as expected[17, 18]. 
This has prompted two different analyses using dispersive methods. They both restrict 
themselves to 7T7r-rescattering but differ in how the subtraction constants are determined 
and in how the dispersion theory is used. Ref. [19] used the Khuri-Treiman equations 
and fixed the subtraction constants by comparing with the one-loop expression at various 
kinematical points. Ref. [20] generalized the reasoning behind [40] to obtain a series of 
dispersion relations for the Mj(t) defined in (3.21). Their predictions on the decay width 
are in agreement within the quoted uncertainties, [20] finds T = 219±22 eV and [19] obtain 
r = 209 ± 20 eV. Both used Dashen's theorem and the then known value of Lg to fix the 
overall constant. These, however, must be compared with T exp = 295 eV which provides a 
check on Dashen's theorem [39]. On the other hand, as emphasized in [11, 49], they have a 
rather different behaviour in the phasespace distributions. This is shown in Fig. 9. It can 
be seen that the slope for [19] is smaller than the order p 4 results while [20] has a larger 
slope than the one-loop result. The latter also follows from the very simplified dispersive 
analysis performed in [11] shown in Fig. 10(a). The general feature of the result of [20] 
seems to be robust against small changes. One motivation for the present work is to check 
these predictions and see if other effects could be important as well, at order p 4 only about 
half the correction came from the unitarity correction. 

We can now compare the same plot coming from our full NNLO calculation shown in 
Fig. 10(b). The total order p G correction is somewhat larger than observed in [19, 20] but 
the trend is definitely in better agreement with [20]. 




In [20] the fact that the position of the Adler zero did not change much from LO to 
NLO was used to determine the subtraction constants. In Fig. 11 we show a blowup of 
the region around the Adler zero. As can be seen, the position of the zero in the real part 
varies a bit by going from LO to NLO and then almost moves back at NNLO. This is not 
incompatible with the observation of [20], we use a slightly different version of NLO then 
the one in [10] and use different input parameters resulting from the order determination 
of LECs [34]. 

The current algebra prediction for the 
Adler zero 6 is s = \m\, independent of 
t. The NLO effect produces a t depen- 15 
dent shift in the zero for the real part. , 
A fairly large positive shift occurs along 
the line t = u going from LO to NLO „ 
and a smaller negative shift from NLO to I 
NNLO as can be seen in Fig. 7(a). Along ' ^ 
the line s = u, a smaller positive shift ap- 
pears from LO to NLO at s = 1.70 m\ 1 
shifting to s = 1.17 ml at NNLO. The . 15 
effects on the slopes can be judged from 

Figs. 7(a) and Fig. 11. 2 o.o2 0.022 0.024 0.026 0.028 0.03 0.032 0.034 

s [GeV 2 ] 

6.3 Dalitz Plot Parameters Figure 11: M(s,t,u) along the line with s = u 

In general we use the data averages of the concentrated along the Adler zero, 
particle data group (PDG)[13]. For the 

distributions in the Dalitz plot the PDG presents no averages. This distribution is usually 
6 We take here the zero of the real part of the amplitude as the Adler zero. 
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described in terms of the variables 



y = —-i 

Qr/ 

Qr) = m v — 2m n + - m^o (6.1) 

for the charged decay. T l is the kinetic energy of pion ir 1 in the final state. The standard 
parameterization of the Dalitz plot is (up to third order) 

|M| 2 = A 2 (1 + ay + by 2 + dx 2 + fy 3 + gx 2 y + •••). (6.2) 

Odd terms in x are forbidden by charge conjugation, all experimental results find them 
compatible with zero and since all most precise experiments have presented fits with the 
odd terms set to zero we use those. / has only been measured by KLOE [50] and no 
experiment has attempted to determine g. Earlier experiments that only determined a and 
b are not included. The results are shown in Tab. 2. There are discrepancies among data 



Exp. 


a 


b 


d 


KLOE [50] 
Crystal Barrel [51] 
Layter et al. [52] 
Gormley et al. [53] 


-1.090 ±0.0051^ 
-1.22 ±0.07 
-1.08 ±0.014 
-1.17 ±0.02 


0.124 ±0.006 ±0.010 
0.22 ±0.11 
0.034 ±0.027 
0.21 ±0.03 


0.057 ± 0.006±g;gS| 
0.06 ± 0.04 (input) 
0.046 ± 0.031 
0.06 ± 0.04 



Table 2: Measurements of the Dalitz plot distributions in rj — > ir + n 7r°. The parameters are 
defined in Eq. (6.2). The KLOE result [50] for / is / = 0.14 ± 0.01 ± 0.02. None of the others 
determined /. The Crystal Barrel used d as input, but remarked that a and b varied very little 
within the range of d used. 

which are hard for us to discuss since correlations are important. There is however a clear 
discrepancy between KLOE[50] and Gormley et al.[53] for a, KLOE[50], Gormley et al.[53] 
and Layter et al.[52] have three rather incompatible numbers for b. The results for d are 
all compatible. 

Similarly the neutral decay is parameterized by 

+ + , = §£(Mz^)\ (6 .3) 

Ei is the energy of the ith pion in the final state. The parameter a has been measured 
by several experiments. There are two recent high precision measurements but they are in 
disagreement, however, KLOE published a new analysis very recently. 

Let us now discuss how to extract the Dalitz parameters in chiral perturbation theory. 
First for the charged decay. The Dalitz plot variables x and y are related to the kinetic 
energy of the outgoing particles. In terms of Mandelstam parameters they are given as, 



((m v - ra^o f - s) - 1 



(s - s) . 



(6.4) 



Exp. 


a 


KLOE [54] 
KLOE (prel)[55] 
Crystal Ball [56] 
WASA/CELSIUS [57] 
Crystal Barrel [58] 
GAMS2000 [59] 
SND [60] 


-0.027 ± 0.004lg;g^ 
-0.014 ±0.005 ±0.004 

-0.031 ± 0.004 
-0.026 ±0.010 ±0.010 
-0.052 ±0.017 ±0.010 

-0.022 ± 0.023 
-0.010 ±0.021 ±0.010 



Table 3 

n 



7r<W 



Measurements of the Dalitz plot distribution in 
The parameter a is defined in Eq. (6.3). 



The first equality is valid in general, the second only in the isospin limit. Our amplitudes 
are in the isospin limit, with a common pion mass everywhere. The physical value gives 
Q v = m v — 2m n + — m n o = 0.13318. For the pion mass we used, 3m 2 = 2m 2 + + m 2 , we 
get Q v = 0.13313 MeV. So the value of is fine. In the isospin limit, x = y = and 
s = t = u = s coincide. In the physical case there is a small difference, s = t = u = sq 
corresponds to y = —0.052. 

We evaluate the decay ampli- 
tude M(s, t, u) in the s-t plane over 
the physical region. Since only 
two variables s and t are indepen- 
dent, the relation s + t + u = m 2 + 
2m 2 . + m 2 Q is used to eliminate 
the third Mandelstam variable u 
in the amplitude. We then fit (6.2) 
to the amplitude |M| 2 = [M^ 2 ) + 
MW + M( 6 )| 2 with as error A = 
ReM( 6 )Re(M( 2 ) + + ) + 
ImM( 6 )lm(M( 2 ) + M< 4 ' + M^). 

This is half of the NNLO contribution. The fitting error quoted in the result (6.5) is from 
this error, not from errors in the input values of the C\ and other LECs. For y we used 
the second expression in (6.4). The fits have been performed to |M( 2 )| 2 , |M^ 2 ) + M^ 4 )| 2 
and |M( 2 ) + M( 4 ) + M( 6 )| 2 labeled LO, NLO and NNLO respectively. In NNLO we have 
studied the effect of setting the C\ = and setting L[ = C\ = and in NLO for L[ = 0. 
In addition, we have checked how much changing y from the second to the first expression 
in (6.4) changes the results, labeled NNLOp as well as including the g term and the terms 
with x 4 ,x 2 y 2 ,y 4 , labeled NNLOq. The results are given in Tab. 4. The changes in going 
from NLO to NNLO are rather modest, the main change is the overall normalization ^4q. 
The linear slope a is lowered somewhat but not enough to reach the most recent experi- 
mental value. The same comment applies to the quadratic slope b. d is in agreement with 
the experimental values. / is always much smaller than the measured value of KLOE. The 
results from the dispersive calculations using Khuri-Treiman[19] 7 and the simplified anal- 
ysis of [11] with two different boundary conditions are also shown in the table. [21] used 
these parameters as input and we have therefore not shown their result. As final result we 
take the NNLO result of Tab. 4 with the MINUIT errors with inputs as above. 

Al = 538 ± 18 , 
a = -1.271 ±0.075, 
b = 0.394 ±0.102, 



7 Note that there exists an analysis [61] using the method of [19] fitting the preliminary KLOE results. 
However, their final normalization depends on the allowed ranges of parameters given by [19]. This is why 
we do not quote their numbers below. 
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KQK 
OOO 


1 OK 1 

— l.zo / 


u.oy / 


U.U/D 


U.UU4 


NNLO [p = 0.6 GeV) 


543 


— 1.300 


0.415 


0.055 


0.024 


1N1NL/U (^/i — U.y \ ) 


D4o 


i i/i i 


U.o ( 4 


U.U04 


n not; 


NNLO (C[ = 0) 


465 


-1.297 


0.404 


0.058 


0.032 


NNLO [L\ = CT = Q) 


251 


-1.241 


0.424 


0.050 


0.007 


dispersive [19] 




-1.33 


0.2b 


0.10 




tree dispersive[ll] 




-1.10 


0.33 


0.001 




absolute dispersive[ll] 




-1.21 


0.33 


0.04 





Table 4: Theoretical estimate of the Dalitz plot distributions in 77 — > 7r + 7r 7r°. The parameters 
are defined in Eq. (6.2). The line labeled NNLO is our central result. 



d = 0.055 ±0.057, 

/ = 0.025 ±0.160. (6.5) 

The fitting results for the neutral decay are shown in Tab. 5. Again we have fitted 
LO, NLO, NNLO and removed the contributions from C\ and L\ to see their effects. The 
one labeled NNLOq is the same fit as NNLO but with the next two terms in the expansion 
around s = t = u = sq included in the fit. As in the case of the charged decay, one sees that 
going from NLO to NNLO does not seem to change much in the Dalitz plot distributions 
but changes mainly the overall normalization. It should be mentioned that because of the 
three terms in (3.20) large cancellations happen in the amplitude for 77 — > 3-7T . This, as 
well as an inequality between the slope parameters is derived in App. A. The inequality is 
always satisfied by our results and the equality which results under additional assumptions 
is reasonably well satisfied. Assuming the error on the NNLO result to be one half the 
NNLO contribution as we did for the charged case leads to 

A 2 = 4790 ± lbO 
a = 0.013 ±0.032. (6.6) 

The results from the dispersive calculations using Khuri-Treiman[19] and the simplified 
analysis of [11] with two different boundary conditions are also shown in the table. No 
theoretical prediction comes near the experimental results of KLOE[54] and Crystal Ball[56] 
except for the result 8 of [21]. 

8 That reference combines a Chiral Lagrangian including the rj and a model for unitary resummations. 
Their good agreement with data is remarkable but should also be understood using more controlled ap- 
proximations. 





-J' 2 


a 


LO 


1090 


0.000 


NLO 


2810 


0.013 




21UU 


U.UlO 


IN IN L(J 


4790 


0.013 


JNJNLOq 


4790 


0.014 


TVTTVTT C\ ( C^T n\ 

INiNL/U ■ — UJ 


a i /in 
414U 


a ni i 
U.U11 


NNLO {Lr = C[ = 0) 


2220 


0.016 


dispersive [19] 




-(0.007—0.014) 


tree dispersive [11] 




-0.0065 


absolute dispersive [11] 




-0.007 


[21] 




-0.031 



Table 5: Theoretical estimates of the Dalitz plot distribution in r\ — > 7r 7r°7r°. The parameter a is 
defined in Eq. (6.3). The line labeled NNLO is the main result. 

It is clear from these results that further study on possible variations of input param- 
eters of ChPT is needed. One puzzling observation is however that the main effect in the 
dispersive calculations is S'-wave rescattering and the input parameters used here give a 
very good prediction for the scattering length a[][62]. 

6.4 The ratio r and decay rates 

We can now proceed to calculate the various decay rates from our amplitude. Since the 
evaluation of the NNLO terms is very slow we have used the following procedure. The actual 
shape of the allowed s,t,u region in the physical charged decay is somewhat different from 
the one allowed with an equal pion mass. We therefore perform the integration over the 
physically allowed values of s, t, u for the neutral and charged case. 

The matrix-element \M\ 2 is dealt with similar to the previous subsection. First we 
fit |M| 2 with an expansion in x and y up to fourth order in x and y for the charged and 
neutral case. The values for x and y are calculated using the values for s, t, u with the 
formulas of the isospin limit with the pion mass as defined above. 

Doing that leads to the decay widths for the decays where we have indicated the various 
orders 

LO, 
NLO , 
NNLO, 
NNLO CI = , 

LO, 
NLO, 
NNLO, 

NNLO(7[ = 0. (6.7) 



nv- 


> ir ir~ir ) = sin 


e 


• 0.572 MeV 




sin 2 


e 


• 1.59 MeV 




sin 2 


e 


• 2.68 MeV 




sin 2 


e 


• 2.33 MeV 


r( v - 


vrVvr ) = sin 2 


e 


• 0.884 MeV 




sin 2 


e 


• 2.31 MeV 




sin 2 


e 


• 3.94 MeV 




sin 2 


e 


• 3.40 MeV 



The numbers in (6.7) can be used to calculate the ratio of decay rates 



- r(r? vr+Tr-vr ) 1 J 

To lowest order and neglecting the differences in phase space this ratio is expected to be 
exactly 1.5. The correct treatment of phasespace and the slightly different pion mass used 
lead to 

r LO = 1.54. (6.9) 
At higher orders a significant change is found with 

rNLO = 1-46. 

?"NNLO = 1-47. 

?~nnloc7=o = 1-46. (6.10) 

The small changes from NLO to NNLO is because the higher order corrections mainly 
change the overall size of the amplitude, not its shape. 

This should be compared with the experimental result [13] of 

r = 1.49 ± 0.06 our average . 

r = 1.43 ±0.04 our fit, (6.11) 

The different results are from the direct measurements or from the global fit to eta decays. 
Our results agree excellently with each value. 

6.5 Discussion and the values of R and Q 

One of the prime reasons to study the hadronic eta decay to pions is the direct determina- 
tion of the double quark mass ratio, Q 2 defined as [39] 

Q 2 = 4^4 (6-12) 

m 2 d -ml 

The reason for considering this quantity is that it is to first order independent of a shift in 
the quark masses of the form 

m u -»■ m u + am d m s , 
m d -»■ m d + am s m u , 

m s — > m s + am u md ■ (6.13) 

Within the framework of standard ChPT such a change is unobservable, it can always 
be compensated by a change in the values of the LECs up to order p 8 effects. This was 
observed to order p 4 in [63] but the generalization to higher orders is obviously true. In 
our case, we are implicitly using assumptions on the values of the constants such that this 
shift is fixed as discussed in [34] . That is one of the reasons we pulled out an overall factor 



of sin(e) rather than Q 2 out of the eta decay amplitude. To first order in isospin breaking 
we have 

• f\ V3m d -m u 1 

sin(e) = e = — — = — — , (6.14) 



m 



4 R 



4 m s 

where the last equality defines the ratio R. 

Since the process is 77 — > 37r is strongly protected from the electromagnetic interactions 
due to the chiral symmetry, we expect to obtain the value ma — m u rather well from this 
process. So let us see what our results imply. Using the experimental value[13] 



F(r, 



7T + 7T 7T°) 



295 ± 17 eV 



(6.15) 



and (6.7) and (6.14) we obtain the values for R quoted in Tab. 6) in the line labeled R(rf). 
For comparison we also quote the values obtained at LO using 





LO 


NLO 


NNLO 


NNLO (C[ = 0) 


R( V ) 


19.1 


31.8 


42.2 


38.7 


R (Dashen) 


44 


44 


37 




R (Dashen-violation) 


36 


37 


32 




Q (v) 


15.6 


20.1 


23.2 


22.2 


Q (Dashen) 


24 


24 


22 




Q (Dashen-violation) 


22 


22 


20 





Table 6: The isospin breaking quantities R are evaluated at p 2 , p 4 and p 6 . The Q values are given 
with Q 2 w 12.7 R. Sec text for details. 



m \-o + m \-+ — 2m 2 o 
R lo = f 5 K+ 2 / (6-16) 
2 (m 2 K0 - m 2 K+ ) 

where the QCD part of the masses should be included only. The electromagnetic part of 
the K° and 7T° mass is taken to be negligible and we use 



m 



K+e 



x D (m^ - m n0 ) . (6.17) 



We take X£> = 1 in the line labeled Dashen[12] and use the estimate xp = 1.84 [64] for 
the line labeled Dashen-violation. The columns labeled NLO and NNLO quote the results 
from [34] using the same input. These used m s /rh = 24 as input but changing it to 26 
changes 9 the value of R by about half a unit. Note that the limit R < 44 [39] is satisfied 
by all the estimates in Tab. 6. 

It is harder to draw conclusions on the value of Q. Its value can be obtained from 



Q 2 = — f--Rto(l + NNLO) , (6.18) 

where again only the QCD part of the masses should be included and m^m 2 are the 
masses in the isospin limit. The NNLO corrections are known [34] and not negligible. 



3 This can be seen from using the number for m u /md and m s /rh from Fig. 3(b) in [34]. 



However, they do depend at present on the input value of m s /rh used in the chiral fits, 
this can be seen in Fig. 3(a) of [34]. Varying m s /m from 24 to 26 changes Q by about one 
unit. The analysis of [39] relied on the NNLO correction being small. 
The relation between Q 2 and R can be written as 



The analysis of [39] leads to a factor of about 12.7 using m s /rh = 24.4, the input value for 
m s /m = 24 to about 12.5 and for m s /m = 26 to 13.5. After taking the square root, this 
range corresponds to an uncertainty on Q of about unit. For completeness we have listed 
in Tab. 6 the values of Q derived from R with the factor equal to 12.7. 

7. Conclusions 

In this paper we have performed a NNLO calculation in standard ChPT of the amplitude 
for 77 — > 3vr. This calculation was performed to first order in isospin breaking and we have 
pulled out the overall factor in the form of sin(e), the lowest order TT -r] mixing angle. The 
remainder of the amplitude is then dealt with in the isospin limit. How we have dealt with 
the pion mass is described in Sect. 6.1 and with the Dalitz plot distributions and decay 
rates in Sects. 6.3 and 6.4. 

We find a reasonable enhancement over the NLO result of [10] which is somewhat 
larger than the estimates from the dispersive analyses [19, 20]. The shape of the amplitude 
is in better agreement with [20] from comparing the published plots along the line s = u. 
We have also commented on the position of the Adler zero which was used in [20] to 
determine their subtraction constants. The NNLO result for the Dalitz plot distributions 
moves somewhat in the direction of the experimental results compared to the NLO result 
but is not in good agreement. The same is also true for the slope parameter a in the 
neutral decay. We always obtain a positive value while experiment is consistently negative. 
The amplitude for this decay has large cancellations and that makes a a difficult parameter 
to predict, however all methods which include unitarity resummations and have published 
values [11, 19, 21] get a negative value. This we find puzzling, since the main effect is tttt 
S-wave rescattering and our input values give a good value for the scattering length afj [62] . 
We obtain very good agreement with the ratio r. 

Since we find a somewhat larger enhancement of the decay rate then [20] we also find 
a somewhat larger value for the isospin breaking quantities R and Q, which means that 
md — m u is somewhat smaller than obtained in [39]. We also find values that are not in 
agreement with the NNLO order fit to the meson masses [34]. 

The influence of changes in the input values is under study, a first impression can be 
had from looking at the results for the C\ = and the different choices of [i. A more 
detailed analysis is planned. 
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A. A discussion on Dalitz plot parameters and the sign of a 

The goal of this appendix is to point a few simple observations about relations between the 
Dalitzplot parameters. We start by parameterinzing the amplitude for the charged decay 
rate as 

M(s, t,u) = A(l + a(s - s ) + b(s - s ) 2 + d(u — t) 2 -\ ) 

= A (1 + ay + by 2 + dx 2 H ) (A.l) 

By computing \M(s,t,u)\ 2 and using (6.4) this leads to the relations 

a = 2Re(a) = -2i?,,Re(a) , 
b = \a\ 2 + 2Re(fe) = R 2 (ja| 2 + 2Re(6)) , 

d = 2Re(d) = 6R 2 Re{d) . (A.2) 

Here we defined R rj = (2m v Q v )/3. We can now use relation (3.20) and s + t + u = 3so- 
We obtain 

M(s,t,u) = A(l + az^ ) 

= A (3 + (b + 3d) ((s - s ) 2 + (t- s ) 2 + (u- s ) 2 ) + • • •) (A.3) 

Using the definition of z these two can be related and give 

A = 3A, (A.4) 

as well 

a = Re(a) = -R 2 Re(b + 3ctj = -(d + b- R 2 \d\ 2 ) . (A.5) 
We thus obtain the relation 

a <\(d + b--A . (A.6) 

Under the assumption that Im(a) = 0, this turns into an equality. If one looks at the 
numbers in Tabs. 2 and 4 we see that there is a very strong cancellation on the r.h.s. of 
(A.6). Note that the KLOE results satisfy the relation with equality quite well and the 
theory results satisfy it within 30% or so. The inequality is satisfied in all cases. The 
underlying reason why it is difficult to get a negative a seems to be that the value obtained 
for b is too large compared to experiment. 

For our NNLO result we can perform the fit also to the amplitude directly and we 
obtain 

A = -22.7 -i 4.38, 
a = -0.631 - i0.183, 
b = -0.017 + i 0.025, 

d = 0.040 - i0.023. (A.7) 



We see that the relation (A. 2) is satisfied. In general, a is sizable but 6, c, and the higher 
order in x, y generalizations are small. 



B. The order p 4 expression 

In this appendix we quote the precise way in which we have defined the order p 4 contribu- 
tion. It agrees with the result of [10]. We use the notation 

P V tv = m 2 - m 2 . (B.l) 

This appears after Am 2 and Am 2 , of (3.11) have been rewritten in the physical masses. 
The functions are defined in (4.2). 

M(5 4) (s) = +L r & (32m 2 K o s - 32 s - 256/3 m 2 m 2 K o + 256/3 m* ) 
+L r 7 (96m 2 K o s - 96m 2 s - 640/3 m 2 m\ a + 640/3 m^,) 
+L r 5 (-8m^ + 64/9 m^o m 2 K0 + 32/9 m 4 ^j 

+L r 3 ( - 16/9 s 2 + 16/3 m 2 K o s - 64/27 m 4 K o + 32/3 m 2 s - 256/27 m 2 m 2 K o - 256/27 m^) 

+ 16^2 ( 1 / 9g2 + 1 / 6m x<> s - 2/27 m^ + 1/3 m^, s - 8/27 m 2 - 8/27 m^ ) 
+A(m 2 )P„- 1 (4/3m 2 oS - 16/9 m^) + ^(m 2 ) (- 3/2 s + 4/9 m^ + 20/9 m 2 ) 
+A(m 2 Y o)^~ 1 (-4/3m 2 oS + 16/9m^ ) +^(m^ ) ( - 1/2 s + 2/9 m^ - 2/9 m 2 ) 
+A(m 2 ) (l/2s-2/3m 2 ) 

+5(m 2 , m 2 o , s) ( - 2/3 s 2 - 4/9 m 2 K o s + 5/3 m 2 s + 2/9 m 2 m^o - 2/3 m£o ) 
+£(m 2 ,m 2 , s) ( - 2/9 m 2 s - 4/27 m 2 m^ + 4/9m* ) 
+B(m^o,m^ ,s) (l/2s 2 - 2/3 m^ s + 4/9m^ - 1/3 m 2 s) 
+B{m 2 K0 , m 2 K0 , 0) (3/2 m^o s - 2/3 m^ - 1/2 m 2 „ s - 2/3 



2 2 \ 



+B(m 2 , m 2 , s) (2/9 m 2 m 2 R0 - 2/9 m£o) 

+C(m 2 K o,m 2 K o,m 2 KO ,s) ( - l/2m^ s 2 + 2/3m^ s + l/2m 2 s 2 - 2/3m 2 m\o s 



(B.2) 



Mf>{t) = J B(m 2 „,m 2 „,t) (- 1/12 i + 1/3 m 2 „) + B(m 2 K0 , m 2 K0 ,t) (- 1/24 1 + l/6m^ ) . 

(B.3) 



M?\t) = Ll (4/3t 2 ) + T ^(-l/m 2 ) 

+B(m 2 , m 2 , t) ( - 1/4 t 2 + 1/3 m 2 ^ i + 1/2 m 2 t - 2/3 m 2 m 2 ^ 

+S(m 2 , m 2 , t) (l/6 m 2 1 - 2/9 m 2 m^ ) + B{m 2 K0 , m^o , *) (3/8 1 2 - m 2 R0 1 + 2/3 m^ ) . 

(B.4) 

C. The order p 6 LECs dependent part 

In this appendix we give the amplitude dependence on the order p 6 LECs C\ in the form 
of Mf (t) as defined (4.4). 

M c (s) = +C*£ 3 (256 m^o s - 512/9 m^o - 128 m 2 m^ s - 1280/9 m 2 ,, m% - 128 m* s 
-5120/9 m^o m^> + 768m£o) 

+CJ 2 (l28m^ s - 512/9 m^o - 64m 2 m^o s - 1280/9 m 2 
-64 s - 1280/9 m^o m 2 R0 + 1024/3 m^ ) 
+C5i (l28m^ s - 512/9 m^o - 64t77 2 m ^ s _ 256/3 777 2 m^ 
-64 s - 1024/3m^o m 2 K o + 4352/9 



+C r 29 (128/3777^0 m 2 K0 s + 64/3 m£o s - 256/9 777^ m 2 R0 - 512/9 m£o 
+C 28 (256/3 m^o s - 512/3 777 2 m^o s - 512/3 m 2 ,) m^ + 256/3 m^ 
+1024/3 777^, m 2 K0 - 512/3 m£o) 
+C r 27 (l28/3 777^o s - 256/3 777 2 o m 2 K0 s - 512/9 m 2 777^0 + 128/3 m£o s 
+1792/9777^0 m 2 K0 - 1280/9 777^0) 

+C r 2& ( - 64/3t77 2 777^0 a + 256/9 777 2 o m K0 - 32/3 m^, s + 128/3 m^ ro^o - 256/9 m£o) 
+C£ 5 ( - 32/3777^0 s 2 + 64/9 mj^, s - 32/9t77 2 „ s 2 + 256/9 777 2 o m 2 K0 s - 256/27 777 2 m K0 
+256/9 777^, s - 1024/27 m^o m^ - 1024/27 m£o) 
+C r 2A ( - 160/9 777^0 s 2 - 64/3 m^o s + 512/27 m% + 160/9 m 2 s 2 

+ 128/3 777 2 o 777^0 S + 256/9 777 2 o 777^0 - 64/3 777^0 S - 256/9 777^0 777^0 " 512/27 777^) 
+ CI2 ( - 64/9777^0 S 2 - 256/9 777^0 S + 512/27 777|-o + 64/3777 2 o s 2 

+128/9t77 2 m^o a + 1024/27 777 2 o 777^0 - 448/9 m£o a + 256/27777^0 m^o + 512/27777^0) 
+C r 2l ( - 256/3 777^„ + 64 m 2 K0 + 64/3 m£o) 



+C20 ^128 777,^-0 s — 256/3 m^ — 64m^ "i^o a — 256/3 m^ m 4 K o — 64m 4 a 
-192m^ m^o + 1088/3 m^o) 

+C[ 9 ( 192 777,^-0 a - 256/3 m^ - 96777 2 m|-o a - 128 m 2 "i^-o - 96m^. a 
-320 m^o m 2 K o + 1600/3 m£ ) 

+C[ 8 ( - 320/3 m 4 K0 a + 512/27 m^, + 352/3 m 2 „ m 2 K0 a + 256/3 m 2 o m 4 K0 - 32/3 m£o a 
- 128/9 m^o - 2432/27 m^) 

+C[ 7 ( - 64/3 m 4 K o s + 512/27 m^o + 32m^ m 2 K o a - 256/27 m 2 ^ m^o - 80/3 m* a 
+1472/27 m^o m 2 ^ - 128/3 m^) 

+C[ 6 (448/3 m^ a + 256/9 m^-o - 704/3 m^o m 2 K o a - 2816/9 m 2 m 4 K o + 400/3 m 4 a 
+4288/9 m^o m 2 K0 - 256 m£„) 

+C[ 5 (512/27 m^o + 256/9m^ m^o — 256/9m^. m 2 K o — 512/27 777,^0 ) 

+C[ 4 ( - 64/3 a + 512/27 m^o - 32t77 2 m\ a + 512/27 m 2 m^o + 112/3 a 

+2240/27 m^o m 2 ^ - 896/9 m£„) 

+C[ 3 ( - 128/3 m 2 K0 s 2 + 128m^-o a - 512/27 m^ - 64/3 777 2 a 2 + 320 m^o m^ s 

-1792/9 m 2 m^o + 128 m£ a - 3584/9 m*,, m 2 K0 - 4096/27 m£„) 

+C12 (64/9 m^o a 2 + 512/9 m^-o a + 512/81 m 6 K0 - 64/3 m 2 a 2 - 256/9 m 2 m 2 K o a 

- 1024/9 m 2 m^o + 896/9 m£ a + 1792/27 m^ - 10496/81 m£o) 

+C[i (128/9771^0 a 2 - 128/3 m^o a + 512/27 m 6 K o + 64/9 m 2 a 2 - 320/3 m 2 m^ a 

+256/3 7/120 m 4 K0 - 128/3 a + 1024/9 m£ ^0 + 1024/27 m£o) 

+C[ ( - 32/9 m? K0 s 2 - 64/9 m 4 K0 a + 256/27 m^ + 32/3 m 2 o a 2 - 256/9 m 2 m 2 K0 a 

+256/27 m 2 o m^ - 256/9 m 4 a + 1408/27 m 4 m 2 K o + 128/9 m^ ) 

+C r 9 ( - 160/9 m^o a 2 + 512/27 m^„ + 160/9 m 2 a 2 - 256/9 m 2 m 4 RQ + 128/3 m£ m% 

-896/27 m^ ) 

+C 8 ' ( - 32/9777,^0 a 2 - 64/9777^0 a + 256/27 m^o + 32/9t77 2 m 2 K0 a - 128/27 m 2 m 4 K0 
+320/9 m^„ a - 128/27 m^ m^o - 128/3 m^,) 

+C£ ( - 64/9m|-o a 2 + 64/3t77^ a - 256/27 m^o - 32/9m 2 a 2 + 160/3m 2 m 2 K o a 
-128/3 m 2 m^ + 64/3 m 4 a - 512/9 m 4 m 2 K0 - 512/27 raj ) 



+C r b ( - 32/3 m 2 K o s 2 + 128/9 m^o s + 64/9 m 2 s 2 + 224/9 m 2 m^ a - 896/27 m 2 m^„ 
-64/9 m^o s - 128/27 m^o - 128/27 m£„) 

+CJ ( - 64/9 s 3 + 128/27 m|- s 2 + 128/9 m 4 K0 s - 512/81 m 6 K0 + 256/27 m 2 s 2 
-256/9m 2 m 2 K o s + 128/9 m£o s + 512/27 m* m^ - 1024/81 m£„) 
+C3 ( - 64/9 s 3 + 128/9 m^o s 2 - 128/9 m^ s + 512/81 m 6 K0 + 256/9 m 2 s 2 
-896/9 m 2 m 2 K0 s 

+512/9 m 2 m 4 K0 - 704/9 m^, s + 3328/27 ro^ m 2 K0 + 5632/81 m^ ) 

+C[ (32/9 s 3 - 64/9 m 2 K0 s 2 + 64/9 m 4 K o s - 256/81 m 6 K0 - 128/9 m 2 s 2 

+448/9 m 2 m^o s - 256/9 m 2 m^ + 352/9 m^, s - 1664/27 m£o m^o - 2816/81 m£o) 

(C.l) 

Mf(t) = 0. (C.2) 

M 2 c (i) = C 2 r 5 (32/3 m 2 1 2 ) + C 2 r 4 ( - 32/3 m^o t 2 + 32/3 m 2 1 2 ) + C 2 r 2 ( - 32/3 m 2 R0 i 2 ) 
+C[ 3 (32m^ 1 2 + 16m 2 1 2 ) + C[ 2 (32/3 1 2 ) 
+C[! ( - 32/3 m^t, 1 2 - 16/3 m 2 1 2 ) + C[ ( - 16/3m 2 K0 1 2 ) 
+C£ ( - 32/3 m^, t 2 + 32/3 m 2 1 2 ) + C 8 r ( - 16/3 m 2 ^ i 2 + 8m 2 t 2 ) 
+C 6 r (l6/3m^ 1 2 + 8/3 m 2 i 2 ) + C r 5 (8/3m 2 „ t 2 ) 

+C4 (l6/3i 3 - 32/9 m 2 K0 1 2 - 64/9 m 2 1 2 ) + C£ (l6/3i 3 - 32/3 m 2 K0 1 2 - 64/3 m 2 1 2 ) 
+C[ ( - 8/3 1 3 + 16/3 m 2 K0 1 2 + 32/3 m 2 i 2 ) (C.3) 
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